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ABSTRACT 



The purpose of the investigation is to find how the natural 
frequency of transverse vibration is affected when the center of gravity 
and center of stiffness are not coincident and to determine if it is 
important in calculating the natural frequency of the fundamental mode 
of vibration of ships . 

Two methods, one an approximate analytic method, and the other 
a lumped mass method, are used to determine the strength of the coupling 
of these modes of vibration. The approximate analytic method involves 
the solution of differential equations from an assumed mode shape and 
then the use of Stodola's equation to compute the change in natural 
frequency. The lumped mass method is a numerical procedure requiring a 
considerable amount, of computation. Both methods are in good agreement 
when the ratio of natural transverse to longitudinal frequency is less 
than one-half. 

The approximate analytic method shows that the change in 
transverse frequency of vibration is a function of the product of two " 
factors. One factor is directly proportional to the area times the distance 
of the center of gravity from the center of stiffness squared and inversely 
proportional to the area moment of inertia. The other factor is a com- 
plicated expression, but which is dependent only on the ratio of trans- 
verse to longitudinal frequency as computed by simpler theory in beams, 
neglecting the offset of the center of gravity from the center of stiffness. 
These egressions have been evaluated and plotted against frequency ratio. 

A simple method of computing the change in transverse frequency of vibration 
is provided. 

When the important parameters for calculating the change in 
frequency of vibration are calculated for a merchant ship and submarine, 
it is seen that the coupling of the fundamental modes of vibration in 
these structures is insignificant. 

In the higher modes of vibration, the ratio of transverse to 
longitudinal frequencies approaches unity. It is recommended that the 
coupling in these higher modes be investigated. 

Thesis Supervisors Frank M. Lewis 

Title: Professor of Marine Engineering 
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NOTATION 



A area of a cross section which will support longitudinal stress 

c arbitrary constant 

E modulus of elasticity 

f ratio of transverse to longitudinal frequency computed disregar ding 

the offset of the center of gravity from the center of stiffness 

F force 

g gravity 

I area moment of inertia of a cross section 

K spring constant 

L length of a "beam 

j i distance between masses in the lumped mass method 

M moment 

m mass 

S stress 

u longitudinal displacement 

V shear 

w frequency 

w q ^ longitudinal frequency as computed disregarding the offset of the center 
of gravity from the center of stiffness. 

w ^ transverse frequency as computed disregarding the offset of the 
center of gravity from the center of stiffness. 

w^. transverse frequency with coupled longitudinal motion 

y vertical displacement 

y distance of the center of gravity from the center of stiffness 

© slope 

XL a number determined by the transverse mode of vibration. It equals 

4.73 for the fundamental mode. 

(aA) weight per unit length 
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INTRODUCTION 



The calculation of the natural frequencies of the longitudinal and 
transverse modes of vibration of any structure have long been considered 
independent.^ In a uniform bar, of course, it is true, but for example, 
in a ship's hull where the center of mass and center of stiffness are not 
coincident, it is not necessarily true. A bending of the ship's structure 
causes some motion in the longitudinal direction of the center of gravity. 
This, in turn, will cause some change in the natural frequency of transverse 
vibration. It is the purpose of this thesis to investigate the strength of 
the coupling of longitudinal and transverse modes of vibration when the 
centers of gravity and stiffness are displaced and determine if it has any 
significance in the calculation of natural frequencies of ships. 



PROCEDURE 



Two distinct procedures are used to find the coupling of longitudinal and 
flexural modes of vibration; one method is an approximate analytic method 
and the other is a lumped mass method. The steps to each method are outlined 
in subsections A and B below; the details are in the appendices. Through 
a numerical example the results of the two procedures are then compared 
for agreement. The entire analysis is then interpreted in the light of 
typical values of significant parameters for a submarine and a merchant ship. 

The investigation is restricted to the fundamental, mode of 
vibration. 

In neither the approximate analytic method nor the lumped mass 
method is the deflection of the system due to shear and rotary inertia 
considered. The transverse motion is due to bending only. 

A. Approximate Analytic Method. 

The general procedure in this method is to find a correction to an 
assumed mode shape from the differential equations of a system in which 
bending causes longitudinal motion and longitudinal, forces cause bending. 

The change in mode shape is then used as a measure of the change in natural 
frequency of transverse vibration. The steps in the procedure which are 
carried out in Appendix I are: 

1. Determine the differential, equations of motion of a system in 
which the center of gravity and center of stiffness are not 
coincident. The first equation is a longitud in a l force balance 
equation and the second one is a bending moment equation. 




2 
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y-0 (2) 

6 dt ax* ax J 

(<xA) = weight per unit length 

A = area of cross section 

g = gravity 

u = longitudinal displacement 

y = vertical displacement 

E = Young's modulus 

y = distance of center of gravity from the center of 
stiffness 

I = area moment of inertia of a cross section 

2) Assume a transverse mode shape of a free- free beam in space. 

jr . =[=os X * + Cosh X X ] (3) 

c = arbitrary constant 

XL = a number determined by the mode of vibration. It is 
equal to 4.73 for the fundamental mode of vibration. 

L = length of the structure. 

3) From the differential equation for longitudinal force balance, 
find the longitudinal motion of the center of stiffness pro- 
duced by the mode shape. 

4) This longitudinal motion causes longitudinal forces which 
produce moments about the center of gravity and hence additional 
bending of the system. Determine the additional deflection 
from these forces using a bending moment equation. 

5) With the results of Step 4 recycle Steps 3 and 4. Hence, a 
first and a second approximation to the additional deflection 



can be obtained. 
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6) Use a method similar to that of Stodola to find the trans- 
verse frequency with coupled longitudinal motion. The ratio 
of the square of the transverse frequency with longitudinal 
motion to that without is proportional to the ratio of total 
deflection without longitudinal motion to the total de- 
flection of the beam with longitudinal motion. 



B. Lumped Mass Method. 

The lumped mass procedure is the Prohl-Jtyklestad sequence method 
for calculating the natural frequency of beams of variable cross section. 
This calculation is modified to include the effect of the displacement of 
the center of gravity from the center of stiffness. 

The calculation is based on the fundamental equations of bending 
and strain in beams: 



M_ 

El 




and 



F 



(AE) 



The beam is divided into sections. At each section an inertia loading is 



added: 



in shear equal to 
2 

m w y 

in bending equal to 
m w y (u - y 0) 

and in longitudinal force equal to 
m w (u - y 6) 



m « mass 



w = frequency 

0 = slope 

1 = length 
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The procedure is to calculate the longitudinal force and displacement, trans- 
verse deflection, slope, tending moment and shear at each section along the 
beam in proceeding from the end of the beam to the mid-point. Various fre- 
quencies are tried until the boundary conditions sire satisfied. For the 
fundamental mode, the slope, shear, end longitudinal motion at the mid-point 
are zero. 
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RESULTS 



A. Approximate Analytic Method 

The equation for determining the frequency of transverse vibration 



by the approximate analytic method is: 

2 _ 2 Y o 

W t = ot Y q + Y 1 + Y 2 



w 



where : 



w 2 _ (22 -j* f ? g I 

ot .4 (aA) 



(5) 



Y o = =[a0os W - 1 • (ffi 2 ) ] 



A ~2 

Y i = ? y c Cos 



2 k 

r "2 D.. X (Cos y D, - 1) 
XL/2 [■ ^ 1 



’ 1 “ D 1 X 



Cos 



w 



(Cos XL/2 - 1) 



(Cosh XL/2 - 1)' 



(1 - D 1 X 2 ) Cos XL/2 (1 + D^X 2 ) CoshXL/2 



( 6 ) 



(7) 



2 k , 2 4 

2 D, X 4 D, X 



Y 2 = (f ) y c Cos XL/2 [( 1 fe-y - — 5 - 5 - 5 ) 



1 - X (1 - D x X ) 




) - <-! 



2 4 

D, X 



-JK 



i - D i 2 x (c°s^r 

(Cos XL/2 - 1] 



(l - q x 2 ) 2 Cos XL/2 



(l + q x 2 ) 2 



» (Cosh XL/2 =1) 1 
Cosh XL/2 J 



(8) 



w^ =transverse frequency with coupled longitudinal motion 



w .-transverse frequency as calculated disregarding the offset 
° ' of the center of gravity from the center of stiffness 



Y = total deflection of the assumed mode shape 
o 
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Y 1 = first approximation to the additional deflection 
Yg 3 second approximation to the additional deflection 

When substitutions shown in Appendix A are made, all the terms in brackets 
in equations ( 7 ) and ( 8 ) are, to an approximation, dependent only on fre- 
quency ratio, f o 



w 



f = 



w 



ot 

of 



where 



w Q £ = longitudinal frequency as calculated disregarding the 
offset of the center of gravity from the center of stiffhess. Specifically, 

11 W t ~ "ot 

D . (ilI2) 2 , (ii) 2 

1 '•fn } 'frr 

Similarly, 

i/a . f 2 
2 

Another important result derived in Appendix A is that the frequency ratio, 
f, is independent of mass per unit length. 

These two important results suggest a rearrangement of the 
frequency equation k, to permit easy calculation of the change in frequency 
for any system. 



w t 

(-M = 

'w . ' 
ot 



Y 1 Y 2 

1 + r + r 

o o 



2 
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This equation can be put in a still more usable form. 




( 9 ) 



where g(f) and h(f ), functions which are only dependent on frequency ratio, 
can be obtained from Figures I and II. 



The function for the first approximation to beam deflection is: 



-2 (— ) 

CO s^f^ w 



(Cos % - 1) 



6(f)- 



2 ^ ;rT )cos (f|) ' (i 



1) . (Cos 1L/2 - 

XL ,, . /XL\2< 



1) 



Cos XL/2 T 
Cosh XL/2J 



•(^) 2) Cos r (1 + (° osh wa) 



j~2 Cos XL/2 



1 - 



The function for the second approximation to beam deflection is: 



2(1=) 



h(f) = Cos XL/2 [( f ^ L 4 









;)(“ 



Cos f 



/- /XL\4.2* ' n „tt 

(1 - (^) ) Cos f^ 



0 - ( 



(hk) 

'ft?' 



,XLx2 

aTi 



)’ 



Sim f 2 + (Cos XL/2 - 1) 

(Cos f|) 2 (1 - (~) 2 ) 2 Cos XL/2 



( ( Cosh XL/2 - 1 ) 

(1 + (|£) 2 ) 2 Cosh XL/2 



j - [2 Cos XL/2 



1 



Cos XL/2 I 
Cosh XL/2 J 



These expressions do not have to be evaluated in finding the change in 

Ay 2 

natural frequency. All that has to be evaluated is w , , w , and — with 

Oh ox J- 

g(f) and h(f) from Figures I and II to use the approximate analytic method. 

Above a frequency ratio of one-half the approximate analytic method 
gives poor results. However, it is still desirable to estimate frequency 
changes at higher frequency ratios. An examination of the relative order 
of magnitude of these two functions suggests an alternate method of deter mi ni n g 
frequency changes from data of only one numerical example of the lumped mass 
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FIGURE I 



FUNCTION FOR THE FIRST APiaOXIMATION 
TO THE CHANGE IN ISAM DEFLECTION 



1.7 



g(f) versus x 



1.6 
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ITGUEE II 



JUNCTION FOR TNE SEC 010) APPROXIMATION 
TO THE CHANGE IN BEAM DEFLECTION 



1=7 
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method shown in Figure IV . This alternate method is based on the fact that 






h(f) « |_g(f )J 



To use Figure IV for other beams use the following procedure: 



w. 



1. Determine — from this figure for the desired frequency ratio. 



ot 

2. This number can be used to solve Equation (10) for x. The 
solution is shown as Equation (11). 

3 . Use this value of x in a new frequency equation with 

for the desired application. See Equation (12). 



w 



(— ) = 
v w ^ ' 



ot 



2 2 

1 + a x + a x 
o o 



*-r -] 

w * ^ 

( ZL) ± - 

W ot 1 + a n x 2 + a, 2 X 2 



( 10 ) 

( 11 ) 

( 12 ) 



where 



w. 



= frequency change obtained from Figure IV 



w , 

ot 

a o - ^ ■ -T5 

Ay 2 

~ for the desired application 



= frequency change for the desired application. 

W ot 



B. Results of the Lumped Mass Method. 

The equations for determining the conditions at each mass along 
the beam are shown below. The zero subscript refers to the conditions at the 
previous section^ the 1 subscript refers to the mass at the point being 
evaluated. The actual numerical value of each of these equations becomes the 
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FIGURE III 



LUMPED MASS SYSTEMS 




THREE-MASS SYSTEM 




ELEVEN-MASS SYSTEM 
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subscript zero term when the next point along the beam is evaluated. 
M 1 = M o + v o ^ + y mi w2 ( U 1 " y 9 i) 



V, - V + 



\ v y x 



m i v j e 

9 1 - 9 o + H- + mr 

M i 2 V i 3 

y l = y o + G o X + 2EI + 6 eI 



U 1 " U o “ K 



F 1 = F o + “l W ^ u i ’ y e i) 



where 



M = bending moment 
V = shear force 
F = longitudinal force 
K = spring constant between masses 
l = distance between masses 

One advantage of the lumped mass method is that it is useful in deriving a 
frequency equation. For the three mass systems shown in Figure III, the 
frequency equation is: 



0 + 2 



m w 



(±1 

V-3 1? 



3 El 




1 

3 EIK 



) 



C. Results of a numerical E xa mple. 

In order to check the validity of the results, a numerical example 
was done using both methods. The change in frequency is shown plotted in 
Figure IV as a function of frequency ratio, f . The plot is for a system of 
these properties: 



-!•+- 




H 






v 

o 



I 

o 



CO 



g H 












'//?f 




O 





o 




o 

CT 



O 

v<) 



O 
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A « 1 


in 2 


I = 1/12 


in' 1 


y - 1/4 


in 


a “ * 284 


lbs/in 


ii 

v_n 





In the lumped mass method the system was divided into eleven masses (Figure III). 
For system of fewer masses the agreement between methods was poor. For example. 



at a frequency ratio of one-fourth: 


W ot 


= -974 


analytic method 


w + 






t 

w . 
ot 


= .961 


eleven mass system 


w. 






t 

w , 


= .950 


nine mass system 


ot 












t 

W 4 . 
ot 


= .910 


seven mass system 


For 


all numerical work, w . and w 

’ ot < 




22.4 la El 




»ot= ? 


V T«aJ 




TT 


I EgA 




W o l = L 


V ToaT 


D* Numerical 


values for 


Marine Application 



The approximate analytic method indicates that there are two sig- 
nificant parameters in calculating the change in transverse frequency, namely, 

Ay 2 

frequency ratio, f, and the geometric factor . 

has been evaluated for two ships: 



The geometric factor 
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merchant ship (3) “x ** 

2 

submarine ( 4 ) s» .017 

These numbers, of course, are dependent on the condition of loading and, 
furthermore, they are valid only for the midship section. The geometric factor 
for the merchant represents an upper limit of this number for the "Gopher 
Mariner". Approximate frequency ratios for these ship types are: 
submarine ( 4 ) f .15 

merchant ship f « .25 
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DISCU5SI0N OF RESULTS, CONCUJSIOKS, AMD RECOMMENDATIONS 

A. Accuracy of Methods 

The approximate analytic method is inaccurate above the frequency 
ratio of one-half. The lumped mass method, on the other hand, can be used 
over the entire range of frequency ratio. 

Reference ( 5 ) states that to ensure better than one percent 
accuracy in computing frequencies by the lumped mass method, thirteen 
masses should be used per wave length. If more masses had been used in the 
numerical example, the agreement between methods at the low frequency ratio 
would be better than it already is. For eleven masses there is only a 1.^ 
difference. 

B. Significance in Ship Vibrations. 

When the constants for the submarine and merchant ship stated in 
Part III above are used in Equation (9) in conjunction with Figures I and II, 
the change in transverse frequency is truly insignificant for the fundamental mode. 
For higher modes of vibration, the frequency ratio is larger and the change 
in frequency might not be so insignificant. 

C . General 

Two general methods for calculating frequency change due to an 
offset of the center of gravity from the center of stiffness are given . 

One method is approximate; the other is exact. Also, a method of estimating 
frequency change by a simple calculation from one lumped mass numerical 
example is explained. The coupling of the transverse and longitudinal modes 
of vibration in any structure can be estimated or obtained exactly by the 
methods explained in the text of this thesis. 
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All the work is restricted to the fundamental mode of vibration. 
The lumped mass method, of course, could be used for higher modes . It is 
recommended that the coupling due to the offset of the center of gravity 
from the center of stiffness for higher modes be investigated. 
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APPENDIX A 

DETAILS OF THE APPROXIMATE ANALYTIC METHOD 

The approximate analytic method, outlined below, is divided into 
steps, the numbering of which corresponds to the numbering in the procedure 
Section IIA above. The description which appears in that section is not 
repeated. 

Step 1: Derivation of the Equations of Motion 

With the sign convention shown in Figure A-l an equilibrium 
force balance on a particle of the system yields this differential equation 
of motion. 



(oA) d 2 / A dy -x . „,d 2 u 



(u + g y) + E/, 



at' 



= o 



(i) 



dx 



where : 



(<xA) = weight per unit length 
g = gravity 

y = y 2 ~ y l = distance of center of stiffness from the center of 
gravity 

E = Young's modulus 

u = displacement in the x- direction of the center of stiffness 
from the equilibrium position 

•-«s 



Letting 



3t e 



( ) = ( ) 



the equation can be rewritten 






g 



dg 



(la) 



dx 



A moment balance about the center of gravity of a s m all particle 
can be written using the sign convention of Figure A- II. 
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FIGUIffi A-l 



LONGITUDINAL FORCE BALANCE 



'IN 




+X 



* 




distance to the cenier of bending 



A i» 



distance to the center of gravity 




O o 





dx }!\. 



•Ml 



+x 
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FICUKE A- II 
MOMENT BALANCE 



(S)A 

4 



M 




o o o 



V + IV 






M + dM 



A ( 3 . f flx) 
> 



SIGN CONVENTION 









+x 



* 



+M 
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dx + M - (M +dM) + (S + ^ dx - S ) Ay = 



Let 



dV 



S3*) *JL ax 

g dt 2 



M d v 

El " 72 
dx 



S = E 



du 

dx 



A = area of the cross section 
Equation (2) can be rewritten: 



ctA S y + 



A , 



St £ 



El + EA T 

dx dx~ 



^y = o 



( 2 ) 



(2a) 



Step 2: Assumed Mode Sha pe -J. r*o 

According to the classical theory of vibration, the equation of 
transverse motion for a uniform beam free in space is: 

y , =[cos Xx ♦ ( g^0 g ) C^h Xx] (3) 

L = length of the beam 
x = distance from mid- length 
c = arbitrary constant 

For the fundamental two-noded mode of vibration: 

XL - L.73 = \/ fs L <3a) 

Step 3; Longitudinal Motion 

Equation 3 is substituted into Equation 2 and the resulting 
differential equation solved for the longitudinal motion, u. The arbitrary 
constants are evaluated by assuming no longitudinal motion at mid-length 
and no stress at the ends. 
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X 

II 

o 


o 

II 

p 


x = L/2 


du ^ 

= o 

dx 



The equation for longitudinal motion becomes: 



u = F, Sin - — + F 2 Sin Xx + F 3 Sinch Xx 



1 v^T 



w 



where 



F,- - D, U 2 [ 2 ° lX 2 J[ C0S 1 

1 1 2 L 1 - D-j 2 X 4J L Cos L/2y/D^ - 



F 2 = 



i - d 2 X 



3 1 + D x 4 2 



D =^- 
1 w 2 (aA.) 



( 5 ) 



D 2 = y c X 



_ _ Cos XL/ 2 

D 3 " Cosh XL/2 



(6) 

( 7 ) 



Step 4 : Additional Deflection Caused by Longitudinal Forces 

In this step the objective is to find the additional vertical 
deflection caused by the longitudinal motion and forces. Again referring 
to Figure A-2, the differential equation to be satisfied is: 



dM ds 
dx dx 



Ay = 0 



where 



El 




„ du 
s = E ~ 
dx 



( 8 ) 



y' = additional transverse deflection 
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Equation (8) becomes 



d 3 y' A - d 2 u 

~ " T y 2 

dx- 3 dx 



(9) 



These boundary conditions are applicable in the solution of equation (9)°. 



II 

o 


X = L/2 


o 

II 

515 


x = L/2 


II 

O 


X 

n 

o 


o 

II 

3 


x = 0 



The result of substituting Equation (4) into Equation (9) and 
solving the differential equation is the equation for the additional deflection 
of the beam caused by the coupled longitudinal motion. The equation is; 

P p 

y' = ^ [ F i Cos \/jr + Cos “ X Cos ^ + Constant (10) 

A measure of the total additional deflection in the beam is obtained by 
evaluating Equation (10) at midlength and at the end. 

y 'x - L/2 

XL/2 [- 1 1 



h ” *'x=0 



Ay. 



c Cos 



■^7? COS^ 



(Cos XL/2 - 1) 
2 



! + D x X 



(Cosh XL/2 = 1) "I 
2 Cosh XL/2 J 



(1 - D 1 X' ) Cos XL/2 

( 12 ) 



Step 5s The Second Approximation to Total Transverse Deflection; 

When Equation (10) is then cycled through Steps 3 and 4 and 
evaluated at the mid- length and end point, the result is a second approximation 



to additional deflection of the beam. 
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2 k 
-2 D, X 4 



A - . 4 r — - , « It D 2 X 1 * COa \fir ' 1 

Y 2 ' ( I> (y) c 003 XL/2 [(— i-5-Tf - ~ 2 - j )( ■■ - y - & ) 

1 - X 1 - X Cos 



, D 1 X w L/2 n Sin L/2\^l" 1 (Cos XL/2 - l) 

-Hi) 4 (1-V 2 ) 2 Cos iL / 2 
Td; 1 

. 1 jfi ( C osh XL/2 - 1) ' 

(1 ♦ Dl x 2 ) s #"= c '* ^ 
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Step 6; Calculation of Natural Frequency. 

With an assumed frequency and deflection, the Stodola procedure 
is to calculate the inertia loading and resultant deflection. If the 
assumed frequency equaled the natural frequency, the resultant deflection 



would equal the assumed deflection. Mathematically; 

2 2 
w, y. = w y 
1 "'l o J o 

2 2 y o 

w, = w — 

1 O y. 



(Ik) 

( 15 ) 



y = assumed deflection 
o 

y^ = resultant deflection 

w^ = assumed frequency 

w 1 = natural frequency of the beam 

When the Stodola equation is applied to the approximate analytical 
procedure, the equation for determining the natural frequency of a system 
with a coupling of longitudinal and transverse motion becomes; 

..2 ..2 Y o 

W t " V ot + Y^+ Y 2 

Y q is Equation (3) evaluated at midlength and at the end 
Y 1 is Equation (12) 

Yg is Equation (13) 
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w^ = natural transverse frequency of the system without coupling 
of longitudinal and transverse motion, 
w^ = natural frequency of the system with coupling of longitudinal 
and transverse motion. 



Additional information can be obtained by some algebraic 
manipulation. From Reference (2), the equations for the square of the 
natural longitudinal frequency and transverse frequency of a uniform beam 
are respectively: 



/„ \2 _ /TT \ EgA 

W o£ " L THU 




22. k 2 Egl 

ISA) 



Therefore, 



f 2 _ W ot 



w 



ol 



22, k 1 I 

2 2 A 

TT L 



(16) 



Equation (l6) gives the very important result that the frequency ratio, f, 
is independent of mass per unit length. 

2 

An examination of the expression X yields another important 

2 

result. Equations (3&) and (5) show that X is equal to: 

2 



D 






gEA 






wf (aA) 



w 2 L(aA7 L 2J l ' ,m 



w 



= ( -#) (^?) 2 



If the change in transverse frequency is small, w^. w ^ 
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Therefore, 




( ^•73 



2 

) 



Similarly, it can be shown that: 



These manipulations show that Equations (12) and (13) are made up to two 
factors; one is and the other is a complicated expression which is 
dependent only on frequency ratio. 
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APPENDIX B 

DETAILS OF THE LUMPED MASS SYSTEM 



The equations in the beam sequence calculation are derived from 
the fundamental equation for bending in a beam. 

El — - M 
dx 



In the case of a beam with the center of gravity and center of stiffness 

coincident, this equation and its integrals, using the sign convention 

of Figure B-I, become; 

,2 

El = M + V x 
dx 2 ° 0 



, £ 

El ^ = M x + V ~ + e, 
dx o o 2 1 



El y = 




where 



V q = shear at the boundary or beam end 

M q = moment at the boundary or beam end 

When these quations are used in a step-by-step procedure, the equations 

for moment, shear, slope, and deflection when evaluated at an incremental 

distance to the right of mass 1, shown in Figure B-I, become: 

M = M + V l 
o o 



V = v o + Mi w y x 



m a 

o 



6 = 9 + rr T 

O El 



V f 

o 

2EI 

.2 



Mf V i J 

y = y 0 + V + A" + sir 
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FIGURE B-I 

SIGN CONVENTION FOR THE LUMPED MASS SYSTEM 
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where 



y Q = deflection at the end 
0 q = slope at the end 



In the case of where the center of gravity and center of stiffness 
are not coincident, these equations are modified and two additional equations 
are added. 



V = V + 
1 o 



9=6 + 

1 o 



u, = u - 
1 o 



V l + 
o 


y w 2 ( Ul - y 


(i) 


2 

m^ w 


y i 


(2) 


M 1 


V l 2 




t' 


o 

2EI 

M l 2 V £ 3 


(3) 


6 1 + 
o 

F 


o . o 

2 ei sir 


(t) 


o 

K 


- 


(5) 


2 

m^ w 


(«i - y 


(6) 



where 



u = longitudinal displacement at position zero 
o 

u^ = longitudinal displacement at position one 



K = spring constant = 



AE 



F q = longitudinal force at the end of the beam 
F^ = longitudinal force at position 1 

y = distance of the center of gravity from the center of stiffness 
In the case of a beam vibrating free in space, the conditions at the end 
of the beam are further specified. 



F = m w (u 
o o o 



j*o> 

a - 



M o ■ h y - m o w y (u 0 - J V 



V = m w y 
o o J o 
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Equations (1) through (6) can he evaluated at other positions along the beam 
using a table as shown in Appendix C. To start the calculation, let 



y 



o 



= 1 



u = u 
o o 



e = e 

o o 



Then a table is completed for each mass along the system using the results 
of the previous table to start the next. The end of the last table, which 
will be the conditions at the mid- point of the beam, will yield three 
equations: 

9 - *11 * a 21 u o + ‘31 9 o 

V ' *12 * a 22 u o + *32 9 o 

u ' a 13 * *23 u o + *33 9 o 

If the proper frequency was used in the calculations, these equations all 

equal zero. The procedure is to let u and V equal zero, solve for u q and 
0 q and substitute into the 0 - equation. When 0 equals zero, the proper 
frequency was used in the calculations. 
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APPENDIX C 

SAMPLE CALCULATIONS OF LUMPED MASS METHOD AT A FREQUENCY RATIO 
OF .25 WITH ELEVEN MASSES 



A =1 In . 2 
I = 1/12 in . 4 

a = .284 lbs. /in. ^ 



L 

f 

y 

i 



m 




El 

i 2 

2EI 

X 3 

Sei 

1 

K 



8.24 in. 



.25 



1/4 in. 



.824 in. 
.60666 x 10 ^ 
3.695^3 x 10 8 

.03295 x 10" 5 

.01357 x 10" 5 

.00372 x 10" 5 

.00274 x 10“ 5 



2 

(lbs. -sec. /in.) 
(radians/ sec . ) ‘ 

(l/lbs.-in) 

(l/lbs . ) 

( in « / lb s . ) 

(in. /lbs . ) 



assume w = .975 w ^ 

Therefore, 

mw 2 = 2.13120 x 10 5 

F = I.O 656 O x 10 5 u - .26640 x 10 5 0 
0 0 o 

M = .26640 xlO^u - .06660 x 10 5 9 
o 00 

V o = 1.0650 x 10 5 y o 

The five-page tabular calculation to reach the center of the beam is on 



the following pages . 
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The last page of the tabular calculation yields three equations which would 
all equal zero if the right frequency was assumed at the beginning of the 
calculation . 

e = 1.27136 y Q + .23484 u q + 2.14740 e Q 

V = 13.68628 y Q + .88648 u q + 23.55147 6 o 

u = .00219 y Q + *92853 u + .01907 e Q 

let y =1 

J o 

u = 0 

V = 0 

Solve for u and 9 and substitute into the 9 equation. 

00 

u Q = .00958 

9 q =-.58148 

9 = .02493 

When the assumed frequency is .95 w 0 -j-> i^e five -page tabular calculation 
repeated, 

9 = - .01893 

The result is plotted in Figure C-l,and thus the natural, transverse frequency 

with offset center of gravity is: 

w = .961 w . 

ot 
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